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Abstract — An inlet manifold or intake manifold is the part 
of an engine that supplies the fuel/air mixture to 
the cylinders. The intake manifold is essential for the optimal 
performance of an internal combustion engine. The objective of 
present paper is to predict and analyze the flow through intake 
manifold of four cylinder spark ignition engine. One of the 
important factors is air flow inside the intake manifold; the ideal 
intake manifold distributes flow evenly to the piston valves. 
Even distribution is important to optimize the efficiency of the 
engine. Hence the flow phenomenon inside the intake manifold 
should be fully optimized to produce more engine power with 
better combustion and further reduces the emission. Firstly 
structural analysis has been conducted to set the thickness and 
material suitability of intake manifold against bursting pressure. 
Three-dimensional inlet manifold was modeled and numerically 
analyzed by using the commercially available FLUENT 
software to study the pressure, velocity and flow characteristics 
inside the runner. The steady state analysis has been carried out 
for three different cases: All runners open, 1st & 3rd runners’ 
open and 2 nd & 4th runners open. The predicted results of total 
pressure loss and total outlet mass flow were discussed. Inlet 
pipe and plenum connection creates a back step geometry which 
causes more total pressure loss due to flow recirculation in 
conventional model. Tapering the geometry is causing more inlet 
mass flow due to reduction in total pressure loss in the plenum 
chamber. 

Key words — Intake manifold, FLUENT, SI engine, Flow 
analysis 

I. Introduction 

An intake manifold is one of the primary components 
regarding the performance of an internal combustion engine. 
An intake manifold is usually made up of a plenum inlet duct, 
connected to the plenum are runners depending on the number 
of cylinders which leads to the engine cylinder. Intake 
manifolds have to be designed to improve engine 
performance by avoiding the phenomena like inter-cylinder 
robbery of charge, inertia of the flow in the individual branch 
pipes, resonance of the air masses in the pipes and the 
Helmholtz effect. Tuning the intake manifold means the 
intake runners are of proper size and length to produce the 
highest possible pressure in the cylinder when the intake valve 
closes. S.Karthikeyan[l] shows, pressure waves for the intake 
manifold is simulated using ID AVL-Boost software, to study 
the internal air flow characteristic for the 3 -cylinder diesel 
engine during transient conditions. The effect of helical, spiral, 
and helical-spiral combination manifold configuration on air 
motion and turbulence inside the cylinder of a Direct Injection 
(DI) diesel engine motored at 3000 rpm. using the CFD tool 
has been presented by Benny Paul [2]. An investigation of 


mixture preparation in the intake manifold of a Diesel 
converted engine into LPG spark-ignition engine operation 
was explained by M. A. Jemni[3] and he proposed two 
manifold shapes in order to test the adequate design in view of 
flow and air-gas homogenization. 



Fig.l Systems of fuel injection [ J_] : a) Single Point Injection, b) 
Multipoint Injection, c) Direct Injection; 1 - Fuel supply, 2 - 
Air intake, 3 - Throttle, 4 - Intake manifold, 5 - Fuel injector 

(or injectors), 6 - Engine 

Jae-soon Lee [4] performs a study for the optimal design of the 
intake system by varying the factors which can influence the 
volumetric efficiency, such as the volume of the plenum 
chamber, the length of the intake manifold and the pipe length 
between the surge tank and the plenum chamber, the 
in-cylinder fluid flow field characteristics of a single-cylinder 
engine is analyzed by B. Murali Krishna[5] to see the effect of 
intake manifold inclination at equivalent rated engine speed 
using Particle Image Velocimetry (PIV) under various static 
intake valve lift conditions. A.Martmez-Sanz[6] develops a 
new design of a high performance intake manifold through a 
combination of CAD and FEM. D. V. Boikov[7] designs the 
intake passage and shows it has a considerable influence on 
organization of the processes of mixture formation and 
combustion in an engine. 

II. Methodology 

A. Governing equations: 

Governing equations: 

The conservative form of mass, momentum and energy 
conservation equations, using Einstein’s summation 
convention over repeated indices, are then given by 

Continuity equation: 

dp d'fpu,-] 

— -f — — - = 0 
df dx,- 
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Momentum equation: 

dp™.) d (p u i u j) 

&t d xj 


dx { 
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According to the Stokes’s hypothesis which assumes that the 
bulk viscosity can be neglected, the shear-stress tensor for a 
Newtonian fluid is given by: 






Energy equation: 






+■ 


where the viscous dissipation & is defined as: 

du; 


B. Steps involved in FLUENT 
In CFD simulation following steps is follow: 

1 . Simplifying the geometry 

2. Setting up the model 

3. Meshing of the model which includes reduction in 



Fig. 1 Intake manifold with features 
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Fig. 2 Manifold Geometry side view 


geometry complexity 

4. Defining boundary conditions 

5. CFD-Post for results 

III. WORKBENCH MODEL 

Geometry is created in ANSYS workbench. 

Table 1. Dimensions of intake manifold 


Section 

Diameter (mm) 

Intake 

38 

Outlet 1 

30 

Outlet 2 

30 

Outlet 3 

30 

Outlet 4 

30 


Table2. Material properties 


Aluminum properties 

Youngs’ modulus 

7e+010N-m2 

Poisson’s ratio 

0.346 

Density 

2710 kg/m 3 

Coefficient of thermal expansion 

2.36e' 5 /°K 

Yield strength 

9.5e 7 N-m 2 


IV. MESHING & BOUNDARY CONDITIONS 

Meshing: The accuracy of the results depends highly upon 
the mesh quality. Thus the choice of meshing scheme (grid 
pattern) is very important for fluent to provide accurate 
results. For doing simulation of the intake manifold model we 
have to do first meshing ,in this technique the flow domain is 
converted or split into various subdomain primitives like 
hexahedral and tetrahedral. Care must be taken to ensure 
proper continuity of solution across the common interfaces 
between two subdomains, so that the approximate solutions 
inside various portions can be put together to give a complete 
picture of fluid flow in the entire domain. We use the 
tetrahedral mesh for this purpose which imposed on model. 
Fig. 3 shows the mesh of intake manifold. 



100 00 


Fig. 3 Mesh of intake manifold 
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Mesh Inflation Layers: 

If we plot a typical velocity profile in the near- wall region, 
we can see that we have a large change in velocity in the wall 
normal direction and it is important to CFD simulation that we 
capture this gradient correctly. To do this, it is required to use 
inflation layer meshing to accurately capture the boundary 
layer region for any wall-bounded turbulent flows. Providing 
a suitable inflation mesh for the geometry is strongly tied to 
the choice of the turbulence model, and the flow field we are 
interested in capturing. 



3600 

Fig. 5 Inflation outlet 


Table 3. Mesh details 


Nodes 

176396 

Element 

s 

575599 


Boundary conditions 

Boundary conditions are essential to do a simulation. In 
this problem inlet is open to atmosphere and at outlet suction 
pressure will act due to pistons down motion. So inlet is 
chosen as pressure inlet and outlet is chosen as pressure 
outlet. Pressure inlet boundary condition needs total pressure 
at the inlet. So from atmospheric condition total pressure at 
inlet is 0 pa (gauge pressure). The default value of reference 
pressure in the operation condition is given as 101325. So 
Gauge pressure (0) = Absolute pressure (101325) - reference 
pressure (101325). In the Turbulent - specification method 
needs to choose Intensity and length scale. Turbulent 
intensity value is assumed as per standard cfd assumption. 
Turbulent length scale value is assumed as 7% of inlet 
diameter as per standard cfd assumption. Remaining is 
needed to keep it as default for this problem. 

Table 4. Boundary conditions 


Pressure at flow 

inlet 

0 pascal 

Intake 

38 

Pressure at runner 

-101325 


outlet 

pascals 

Hydraulic 

diameter 

60mm 

Turbulence 

intensity 

4.2% 

Pressure at runner 

outlet 

-101325 

pascals 


In the Turbulent - specification method needs to choose 
Intensity and length scale. Turbulent intensity value is 
assumed as per standard cfd assumption. Turbulent length 
scale value is assumed as 7% of inlet diameter as per standard 
cfd assumption. Remaining is needed to keep it as default for 
this problem. 

V. RESULTS AND DISCUSSION 


A. CASE 1: All runners open 


I I 32e+02 
1 26e+02 
1 19e+02 
1 13e + 02 
1 06e+02 
9 96e«01 
9 29e+01 
8 63e+0l 
7 9/e401 
7.31e+01 
6 65e401 
S99e401 
b33e401 
4 67e401 
4 01C401 
335e + 01 
2 69e401 
2 03c + 0t 
■ 1 37e401 

7 12e400 
I 5 23c 01 


Velocity Vectors Colored By Velocity Magnitude (m/s) Feb 26, 2015 

ANSYS Fluent 15.0 (3d, pbns, ske) 

Fig. 6 Velocity vectors (all runners open) 


| 531e + 01 
-5266402 

1 10e + 03 
1 68e403 
2 26c+03 
2 8-10403 
3 42e.03 
4 00e403 
4 58e403 
-5 166403 
-5 74e»03 
6 32e»03 
6 89e + 03 
-7 47ef03 
-8 05C403 
-8 63e»03 
9 21c + 03 
-9 /9e«03 
• -1 04et04 
1 09e404 
I 1 156404 


Contours of Total Pressure (pascal) Feb 26, 20 15 

ANSYS Fluent 15.0 (3d, pbns. ske) 

Fig. 7 Pressure contours (all runners open) 

Fig. 6 & 7 shows the velocity and pressure contours for all 
runners open. It is observed that velocity drops as the flow 
proceeds through the plenum chamber. This is due to sudden 
increase of the area within the plenum. There is a drop in 
velocity at the inlet of the runner 1 compared to other runners, 
due to a sharp bend at the region of runner 1 inlet. Due to the 
stagnation of fluid occurring at the end corners of the plenum 
high pressure regions are created. 
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B. CASE 2: 1st & 3rd runners open: 


When 1 st and 3 ld runners are open, then the other two runners 
i.e 2 nd and 4 th are considered as wall in the named selection so 
that there will be no flow in those runners. 



1 21c+02 
1 15C+02 
I 09e+02 

1 03© . 02 
9 69o+01 
9 09e *0 I 
8 48e .0 1 
7 88e+0l 
7 2/e. 01 
6 6/e. 01 
6 06e.0l 
5. 450+01 
4 850 +01 
4 240+01 
3 64e.01 

3 03c +01 

2 42e+01 
1 82e+01 
1 21e+01 
6 06c+00 

4 40c 05 



Velocity Vectors Colored By Velocity Magnitude (m/s) 


Feb 27, 2015 
ANSYS Fluent 15 0 (3d, pbns, ske> 


rH 

Fig. 8 Velocity vectors (1st & 3 runner open) 



3 45e+02 
-1 89e+02 
-7 23e+02 
-1 26e.03 

1 79e+03 
-2 33e+03 

2 8Ge+03 



8 20e » 03 
-8 74e+03 
-9 27e+03 
-9 80e+03 


1 03e+04 




Contours of Total Pressure (pascal) 


Feb 27, 2015 
ANSYS Fluent 15.0 (3d, pbns, ske) 


Fig. 9 Pressure contours (1st & 3 ld runner open) 

From Fig. 8 & 9 it is observed that when runner 1 st and 3 rd are 
open the velocity distribution within the plenum changes 
drastically with higher velocity occurring at inlet region of 1 st 
and 3 rd runner. The pressure is higher when 1st and 3rd 
runners are open compared to 2 nd and 4 th runner open 
condition in the plenum. The pressure within is lower 
comparatively for the condition when all runners are open. 


CASE 3: 2nd & 4th runners open 


I I 24e+02 
1 18e+02 
1 lie. 02 
1 05e+0? 
9 91e+01 
9 29e+01 
8 67e.01 
8 Obe.OI 
7 43e.01 
6 Ble.OI 
6 19e+01 
5 57e+01 
4 95e+01 
4 33e .01 
3 71e.01 
3 10e+01 
2 48e.01 
1 86e*01 

H I 24e+01 
6 l9e+00 
1 97e-04 



Velocity Vectors Colored By Velocity Magnitude (m/s) Feb 27, 2015 

ANSYS Fluent 15.0 (3d, pbns, ske) 



Contours of Total Pressure (pascal) Feb 27, 2015 

ANSYS Fluent 15.0 (3d, pbns, ske) 

Fig. 1 1 Pressure contours (2 nd & 4 th runner open) 

Fig. 10 & 11 shows the velocity and pressure flow when 2 nd 
and 4 th runners in open condition. It is observed that when 4 th 
runner is open the velocity increases at the inlet region and the 
pressure is low when 2 nd and 4 th runners are open compared to 
1 st and 3 ld runner open condition in the plenum. The pressure 
within is lower comparatively for the condition when all 
runners are open. 



•all open 
*1&3 open 
•2&4 open 


Plenum length (mm) 


Fig. 12 Velocity variations with plenum length 


From Fig. 12 it is observed that a sudden drop in 
velocity occurs for all the three conditions as the flow 
occurs from inlet of the plenum to the runner. 



Plenum length (mm) 


all open 
1&3 open 
2&4 open 


Fig. 13 Pressure variations with plenum length 

From Fig. 13 it is observed that the pressure is higher when 1 st 
and 3 ld runners are open compared to 2 nd and 4 th runner open 
condition in the plenum. The pressure within is lower 
comparatively for the condition when all runners are open. 


Fig. 10 Velocity Vectors (2 nd & 4 th runner open) 
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Fig. 14 Total pressure of different geometries 
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Total pressure loss 


Fig. 15 Total pressure (1&3 open) 
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■ straight 

■ Top & bottom tapper 

■ Back tapper 

■ Both_tapper 


Fig. 16 Total pressure (2&4 open) 

When compare with base model as well as all modified 
models 

- Due to less total pressure loss, inlet mass flow rate of 
modified geometry of all taper is higher than others. Tapering 
is not causing any impact on the recirculation zone reduction. 
So except Runner 1 remaining Runners are improved in 
carrying the mass flow 


Conclusions 
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After analysis, the concluded points are as follows: 


• All the cases show that Runner 4 is working efficiently 

because it’s located away from the pressure loss or 
recirculation region. 

• Total pressure loss for case 1 (all runners open) is 

8338 pa, for case 2(l st and 3 ld runners open) is 5827 
pa, for case 3 (2 nd and 4 th runners open) is 5500 pa. 

• Outlet total pressure for Case 1 (all runners open) is 

runner 1 is -10543 pa, runner 2 is -8415 pa, runner 3 
is -6896 pa, runner 4 is -7499 pa. Case 2(l st and 3 rd 
runners open) is runner 1 is -6554 pa, runner 3 is 
-4807 pa. Case 3 (2 nd and 4 th runners open) is runner 
2 is -6193 pa, runner 4 is -4769 pa. 
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